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Constraining white-dwarf kicks in globular clusters 
III. Cluster Heating 
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ABSTRACT 

Recent observations of white dwarfs in globular clusters indicate that these stars may get 
a velocity kick during their time as giants. This velocity kick could originate naturally 
if the mass loss while on the asymptotic giant branch is slightly asymmetric. The kicks 
may be large enough to dramatically change the radial distribution of young white dwarfs, 
giving them larger energies than other stars in the cluster. As these energetic white dwarfs 
travel through the cluster they can impart their excess energy on the other stars in the 
cluster. This new heat source for globular clusters is expected to be largest during the 
clusters' youth. 
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1 INTRODUCTION 



ISpruitJ (1 19981 ') proposed that white dwarfs can acquire their 
observed rotation rates from mild kicks generated by asym- 
metric and off-centered winds toward t he end of their time 
on th e asymptotic Riant bra nch (AGB) l|Vassiliadis fc Wood! 
ll993l ). lFellhauer et all (|2003l ) invoked these mild kicks to ex- 
plain a putative dearth of white dwa r fs in open cl u sters (e.g. 
IWeidemannlll977l : iKalirai et al.ll2~00ll ). iDavis et all l|200ct ) ob- 
serv ed a possible sig n ature of white-dwarf kicks in NGC 6397, 
and ICalamida eta l. (2007) fo und si milar but weaker hints in 
Omega Centauri. Davis et alJ l|2008h found that young white 
dwarfs are less centrally concentrated than either their progen- 
itors near the top of the main sequence or older white dwarfs 
whose velocity distribution has had a chance to relax. 

Early in the life of a cluster mass loss due to stellar evo- 
lution com petes with the loss of stars due to evaporation from 
the cluster l|Spitzerlll987l ). Without a kick young white dwarfs 
would have a velocity distribution nearly equal to that of their 
more massive progenitors on the main sequence. In this case 
the kinetic energy of these white dwarfs is much less than 
equipartition; therefore, as their velocity distribution relaxes 
they cool the rest of the cluster. If on the other hand white 
dwarfs receive a substant ial kick at birth as observations in- 
dicate (|Davis et al.11200 8). young white dwarfs may heat the 



rest of the stars in the cluster. This letter will examine how 
white dwarf kicks affect the energy balance within a globular 
cluster. 



2 CALCULATIONS 

Clusters of stars can typically be modelled with a lowered 
isothermal profile (or Kin g model) (|Michiel fl963l : lKinel [l966; 
iBinnev fc Tremaine|[l987l ). 
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where e — ^> — and is the gravitational potential. Be- 
cause the distribution function depends only on constants of 
the motion (the energy), it is constant in time as well. 

With time the kinetic energy within the cluster ap- 
proaches equi partition betw een the various stars such that 
miO~i = rnj <j] ijSpitzerll 19871 ). The progenitors of young white 
dwarfs will be the most massive main-sequence stars in a clus- 
ter at the time, so they will typically have <tto < ^cluster, 
where o" c i ustor is the mean velocity dispersion of the cluster. 
As these stars evolve they lose mass. During a globular clus- 
ter's youth, this stellar mass loss dominates the mass loss from 
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Figure 1. The total binding energy of a cluster modelled with a 
lowered isothermal distribution function as a function of the central 
gravitational potential. For small $(0), the total binding energy 
increases linearly with the central potential. For large values of $(0), 
\E\ oc $(0) 2 . 
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where dN/dm is the number of stars per unit mass in the 
cluster (the mass function), M is the total mass of stars in 
the cluster, m is the mean mass of stars in the cluster, r(m) 
is the duration of the main sequence for a star in the cluster, 
t = r(mTo) (the age of t he cluster), and £ = ("Ukick/c) 2 ~ 
0.85. llben fc Renzinil l)l985n give the initial-final mass relation 
in t he numerator of the rightmost expression. 

iRicher et all i|2007l ) estimate the current mass function 
for a region near the half-light radius of NGC 6397 to be 

dN _ a d\nN 

— = Am , — = 1 - a (3) 

dm d In m 

where the slope of the mass function today is given by a — 0.13 
and jtito ~ O.8M (to « 12 Gyr is the current age of the 
cluster) and dlnr/dlnm w -3.75 (|Straniero fc Chieffi|[l99ll ) 
for a cluster of the age and metallicity of NGC 6397. 

Integrating Eq. over the masses of the stars in the 
cluster gives the mean mass of a star in the cluster (assuming 
a < 1 and mxo 3> m m i n , the minimum mass of a star), 



-771TO 



0.37M 



■ 



CkickT 



it 

2 



'O.38M 



+ 0.15 



• 0.13 



Ma 2 2 ft \ m T o 
for the current observations of NGC 6397. 



(4) 



(5) 



the cluster; later the relaxation of the stellar velocity distribu- 
tion to a Maxwellian and the subsequent evaporation of stars 
from the cluster dominates. 



2.1 Kicks 

At middle age the mass of white dwarfs that remain in the 
cluster may exceed a third of the total mass of the cluster; con- 
sequently, if white dwarfs recei ve a kick comparable to the ve- 
locity dispersion of the cluster (|Davis et alj|2008l ; lHevlll2007al . 
I2008T I. the total kinetic energy of the kicks may approach 10- 
20% of the binding energy of the cluster as shown in Fig. [T] 

Specifically, iHevll <|2008) found that an initial distribution 
of white-dwarf progenitors with <jto = 0.5cr and a typical 
kick velocity of at = 1.84<tto could explain the observations 
|Davis et al.ll2008l ). Essentially, young white dwarfs receive a 
kick of the same order as the velocity dispersion of the cluster. 
Furthermore, unlike neutron stars whose kicks nearly always 
cause them to leave the cluster, most of the young white dwarfs 
remain in the cluster to heat it up (only about two percent 
escape within a crossing time). 

These results provide an estimate of the total power in 
white-dwarf kicks, 



2.2 Binaries 

In a binary a fraction of the energy of th e kick is used to 
change the orbital parameters (|Hevlll2007br ). If the masses of 
the primar y and secondary are simila r due to dynamical bi- 
asing (e.g. iMcDonald fc Clarke! [l993n . the total kick to the 
binary is about 70-80% of the kick imparted to a single star, 
so given that the fraction of binaries is small in the cluster 
as a whole and the correction to the kick for binaries is also 
small, the energetics of changing the orbits of binaries will be 
ignored. 

On the other hand, the binaries provide an important en- 
ergy source for the cluster (in fact the only energy source if one 
excludes the kicks). The interaction of a binary with a single 
star can result in an exchange or the dissolution of the binary, 
but generally it results in an increase in the kinetic energy of 
final single star (like a kick) at the expense of the increased 
binding energy of the binary. The velocity increment of the 
singleton is generally random (also like a kick), so these two 
energy sources are similar and possibly comparably important 
to t he evolution o f the cluster. 

ISpitzerl |l987l ) gives an estimate for the power from bina- 
ries of 

0.92 Ar mcr 2 ... 

Ebinarytr ~ ~wTT b ~o~ W 

per relaxation time, 
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where Nb is the number of binaries. Combining these results 
yields 
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where the binary fraction is taken to be around a few percent 
suitable for a typ ical region of the cluster outside the core 
l|Davis et al.ll200g ). The power produced by the binaries is of 
course proportional to the number density of the binaries and 
inversely proportional to the velocity dispersion that sets the 
cross-section per binary. 

Comparing Eq. (JHJ) to (JSJ shows that the pow er sources are 
simila r for the region of NGC 6397 observed by iRicher et al.l 
(|2007h . 

2.3 Evolution 

To look at the relative importance of binaries and kicks early 
in the life of the cluster, some assumptions about the mass 
function of the cluster in the past are needed. First, the mass 
function becomes more and more top heavy with time as the 
low mass stars are lost from the cluster; therefore, it is natural 
to assume that a > 1 in the past and possibly a > 2 near the 
turnoff. In this regime, the derivation of Eq. ([5]) is not valid. 
The result in general is 
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where the v arious slopes a r e eva luated at the turn-off. Using 
the IMF of iKroupa et all (| 19931 ) below one solar mass and 
IScald (1 19861 ) above gives |1 — a\ «s 1.7 compared to 2 — a w 
1.9 currently. Also the value of the /3 only changes with time 
slightly (for tuto ~ lM©, (3 w 3), so the bulk of the increase 
comes from the replacement of the turn-off mass with the mean 
mass in the denominator of the expression. 

Of course as the cluster evolves, the velocity dispersion of 
the cluster should also evolve. Because the mass of the cluster 
was larger in the past, one would expect that the velocity 
dispersion was also larger. On the other hand, the kick velocity 
may also change with the turn-off mass, so it is natural to 
introduce both of these quantities as variables. The power from 
binaries also depends sensitively on the velocity dispersion and 
mean stellar mass. Taking the ratio of the kick power to the 
binary power yields, 
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where the approximation holds for r = 10 7 — 10 9 yr. After 
about a billion years, one would expect the cluster to have 



evolved structurally, ejecting many of the low mass stars. This 
would change the mass function and typically decrease this 
ratio further. The increased velocity dispersion of the cluster 
actually increases the relative importance of kicks early in the 
life of the cluster by increasing the relaxation time (decreasing 
the binary power). On the other hand, the number density of 
binaries was likely to be larger in the past than today simply 
because the number of density of stars was larger then. 



3 CONCLUSIONS 

Over the life of a globular cluster such as NGC 6397, white- 
dwarfs kicks may provide a significant energy source. In a re- 
gion outside the core todayt, kicks provide about about one- 
half of energy input from binaries. Early in the life of the glob- 
ular cluster when stars of several solar masses are leaving the 
main sequence, the white-dwarf kicks may actually dominate 
over binaries as an energy source; consequently, young glob- 
ular clusters such as those in starburst galaxies may a c tually 
differ structurally from their older peers. [ Drukier et ail (|l992l ) 
found that an energy source beyond binaries was required to 
avoid core collapse in M71; perhaps white dwarf kicks could 
explain this discrepancy. Regardless, these calculations indi- 
cate that further study of the effects of white-dwarf kicks on 
the dynamics of globular clusters is warranted. 
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